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■ Abstract: We discuss the end of the first order phase transition and the bound 
state spectrum, both at weak transition and at the crossover. 

1 Introduction 

Two years ago |], we have presented here first results of our group obtained within 
the 3D approach, comparing the properties of the thermal phase transition in the 
SU (2) Higgs model at small coupling (Higgs mass) = 35 GeV with a more 
realistic = 70 GeV.0 The endpoint of the phase transition and the physics near 
to the corresponding temperature T c {m% %t ) is the topic of this talk. 

In the meantime, due to efforts of three groups 0, SI using the 3D approach, 
various aspects of the high temperature electroweak phase transition (as latent heat 
and interface tension) have been explored in this model within this span of Higgs 
^ ■ masses where the character of the transition changes drastically. The accuracy is 
inaccessible to 4D Monte Carlo simulations [Q, although the results are consistent 
with each other where they can be compared || [F]. 

The interest in the properties of this phase transition resulted from the hope to 
work out - within the standard model - a viable mechanism for the generation of 
baryon asymmetry of the universe at the electroweak scale (see Z. Fodor ||). It 
turned out, however, that in the case of the standard model, taking the lower bound 
of the Higgs mass into account, it is at most weakly first order. Together with the 
small amount of CP violation in the standard model this has ruled out the most 
economic scenario of BAU generation. 

Using dimensional reduction || the simplest effective 3D SU(2) Higgs theory has 
the action 

5-8 = 1 d 3 x(±F b ap F b a/3 + (D^+iD^) + ml<j> + cf> + A 3 (0 + 0) 2 ) (1) 



lr Va\k at Buckow Symp. 97 given by E.-M. Ilgenfritz, supported by DFG under grant Mu932/l-4 
2 Mff, together with the 3D gauge coupling, is labelling our lattice data. For relation to the 
zero-temperature Higgs mass, see Ref. H- 
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with dimensionful, renormalisation group invariant couplings g\ and A3 and a running 
mass squared 777,3 (yU 3 ). ^ ^ s related to the corresponding lattice model with the action 

S = Pa (1 - itrl/p) - fa E itr(^^**f«) + E 0£ + Mp* - I) 2 ) (2) 

p x,a x 

and p 2 = |tr($+$ 2 .) through the following relations (clarifying the meaning of M#) 

^ _A 3 /3|_1/ M£ \ 2 $h * _ 2(1 - 20r) 

agV glPa 8v80GeVy (3 G ' P " 6 + a 2 mi " 1 J 

The coupling parameters in the lattice action can be expressed in terms of 4D cou- 
plings and masses. The bare mass squared is related to the renormalised 777K/Z3) 
(we choose /i 3 = g%) through a lattice two-loop calculation |TIJ. Most of the 3D 



numerical investigations have been done for the SU(2) Higgs model. The results 
obtained in the 3D approach indicate the validity of the dimensional reduction near 
to the transition temperature for Higgs masses between 30 and 100 GeV. 

Our simulations have been done on the DFG-sponsored Quadrics QH2 in Biele- 
feld, on a CRAY-T90 and a Quadrics Q4 at HLRZ in Julich. Algorithmic details are 
described in Ref. [Q. In the course of time, we have explored the parameter range 
= 35, 57, 70, 74, 76, 80 GeV, looking for the strength of the transition, checking 
the approach to the continuum limit by varying /?G = 8,12,16ocl/a and comparing 
with perturbation theory. The simulations at the last three parameter values Mo- 



have been devoted exclusively to locate the endpoint of the transition line |Tl] in 
the rriH — T plane. The lattice size usually varied from 20 3 to 80 3 . In order to 
unambiguously locate the critical Higgs mass we had to increase the lattice to 96 3 
simulating at = 74 GeV, the nearest to the critical value. 

It has been essential to make extensive use of the Ferrenberg-Swendsen multihis- 
togram technique. Data obtained at various Mfj and (3h have been subject to global 
analysis by reweighting to other values, eventually extending the hopping parame- 
ter to complex values. The method constructs an optimal estimator of the spectral 
density of states (at fixed L and Pa) occurring in the representation of the partition 
function 

Z(L, p G , p H ) = J dS 1 dS 2 D L {S 1 , S 2 , p G ) exp (L\p H S x - P 2 H S 2 )) (4) 
from the measured double-histograms in the two variables 

5i = 3E^(7,A +a ), ^ = -^-(Ep'-2Ep')- (5) 

x,a * {JsyG x x 

2 Finding mg rf 

We have employed two methods to localise the critical Higgs mass. First we tried an 
interpolation of the discontinuity 

A(rt)/^ 2 = l/8/5 G ^A(p 2 ), (6) 
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based on reweighting simulation data taken at = 70, 74 and 76 GeV, in order 
to see where this turns to zero. This discontinuity is directly related to the latent 
heat. The pseudocritical (3h c has been defined in two ways, by the maximum of 
the p 2 susceptibility and by the minimum of the p 2 Binder cumulant. For each 
Mff, an infinite volume extrapolation must be performed. For this purpose we have 
expressed the physical lattice size by the dimensionless variable Lag 2 = 4L/ (3q. We 
have collected data from different lattice volumes L 3 , from different /3q (= 12 and 16) 
as well as from different definitions of @h c and found them scattered along a unique 
function of this variable (Fig. p. We tried to fit the finite volume scaling behaviour 
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as A((f) + (f)) OD — A(0 + 0)l oc (L a g 2 )~ 2 suggested by the Potts model. Using only 
the largest lattice volumes available, the criterion of vanishing latent heat gives an 
upper bound \^ lt / g 2 < 0.107 for the existence of a first order transition. 

The other method detects the end of the transition line using the finite size 
analysis of the Lee- Yang zeroes. A genuine phase transition is characterised by non- 
analytical behaviour of the infinite volume free energy density. This is caused by 
zeroes of the partition function (in our case as function of (3h extended to complex 
values) clustering along a line nearest to the real axis and pinching in the thermo- 
dynamical limit. If there is a first order phase transition the first few zeroes are 
expected at 
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(7) 



The partition function for complex couplings is obtained by reweighting from mea- 
surements at real couplings. The first zeroes can be well localised using the Newton- 
Raphson algorithm. We fit the imaginary part of the first zero for each available 
length Lag 2 according to Im/?j^ = C(Lagl)~ L/ + R. The scenario suggested is the 
change of the first order transition into an analytic crossover above M^ rtt . A positive 
R signals that the first zero does not approach the real axis in the thermodynamical 
limit. A similar investigation has been performed recently at smaller gauge cou- 

0.102(2). 



pling in Ref. [12]. The fitted R crosses zero (Fig. 
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corresponds to a zero temperature Higgs mass m% lt = 72.4(9) GeV, and the phase 
transition line ends at a temperature of T c = 110(1.5) GeV. This refers to the exper- 
imental top mass. For the purely bosonic Higgs model we get m% %t = 67.0(8) GeV 
and T c = 154.8(2.6) GeV. 

3 Bound states below and above rrfff 1 

To study the ground and excited bound states with given 

jPC 

one has to use cross- 
correlations between operators Oi forming a complete set in that channel. They 
describe a spectral decomposition tyf 1 ' = (vac|0 i |\I ,(n )) with |\I> (n )) being the (zero 
momentum) energy eigenstates. According to the transfer matrix formalism, the 
connected correlation matrix 



C^t) = (0i(*)Oi(O)) = £ *l n) *J n) *e- m »* (8) 

71=1 



at time separation t gives the masses and wave functions. 

Practically, only a truncated set of operators Oi(i = 1, . . . , N) is used assuming 
that this allows to find the lowest states from the eigenvalue problem for Cy(t). 
However, truncation errors are not small. Solving instead the generalised eigenvalue 
problem 

E^i(^i n) = AM(Mo)E^oK n) (9) 
j j 

these errors can be kept minimal (t > to, to = 0, 1) [13]. 

The wave function of state n at fixed small time t > to is found to be (t) = 
(vac | Oi e _Ht |\I ,< ' n ' ) ). The components in the operator basis characterise the coupling 
of Oi to the (ground or excited) bound state n in the J PC channel. The masses of 
these states are obtained fitting the diagonal elements ^ n \t) 

^\ t ) = j:^ n> c l ,mf (io) 

to a cosh form with t in some plateau region. In contrast to a blocking procedure 



used in jnj, our base is built by a few types of gauge invariant operators Oi, properly 
chosen with respect to lattice symmetry and quantum numbers, having different well- 
defined transverse extensions. In the Higgs channel (0 ++ ) we use the operators p 2 x and 
<3x,n(l) = \^(&tUx^ ■ ■ ■ Ux+Q-i^nQx+ip) as well as quadratic Wilson loops of size I x 
/, in the W-channel (1~~ ) the operators V^(l)) = \ti(T b §+U x ^ . . . C4+(j-i) MlAI $x+Z/<) 
and in the 2 ++ channel S X ^{1) — S Xjl/ (l) where / expresses the operator extension in 
lattice units. In our procedure, the operator $( n ) projecting maximally onto state 
n is a optimised superposition = Yl,i o!^Oi with coefficients provided by the 
(properly normalised) solutions of the generalised eigenvalue equation. The 
coefficients give direct access to the spatial extension of the states under study. This 
construction, at different (5q, should reveal the same wave function in position space 
as function of lag 2 . 
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As examples we present in Figs. [|,^| squared wave functions in the ++ channel 
vs. lag\ near the end of the phase transition at = 70 GeV (from a 30 3 lattice). 
In the symmetric phase the second excitation is a pure H^-ball state in agreement 
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Figure 3: ++ ground state, left: Wilson Figure 4: Same as Fig. 3 for second/first 
loop projection, right: S X:fl (l) projection, excited state 
upper /lower plot: symmetric/Higgs phase 

with jnj] (the first excitation is a Higgs state, too). In the broken phase, the first 
excited Higgs state still contains some admixture of pure gauge matter which vanishes 
only deeply inside this phase. This seems to be a precursor of the near end of the 
transition. More details, statistics and results for other quantum numbers will be 
discussed elsewhere £[5] . 

Since Wilson loops have different projection properties in both phases we have 
studied the flow of the mass spectrum with (3 H or m\{g1)/g\ at = 100 GeV 
(A 3 /(?3 = 0.15625) safely above the critical Higgs mass, too, passing there an ana- 
lytic crossover. In Fig. we present the ++ spectrum in the vicinity and far from 
the crossover. The wave functions are similar to those discussed before. Here the 
spectrum on the "symmetric" side shows decoupling of Higgs and PF-ball states, too, 
the lowest VF-ball mass is roughly independent on (3h- Passing the crossover (from 
"symmetric" to the "broken" side) the first Higgs excitation contains also excited 
glue. These contributions vanish at lower temperature as indicated by the right fig- 
ure. The minimum of the ground state mass remains finite, therefore the correlation 
length fits into the lattice. 

The phenomenological interest is now concentrated on extensions of the stan- 
dard model, with MSSM being the most promising candidate. Concerning non- 
perturbative physics in general, however, the lattice version of the standard Higgs 
model remains interesting as a laboratory for investigating the behaviour of hot gauge 
fields coupled to scalar matter, for the characterisation of possible bound states, for 
the understanding of real time topological transitions and as a cross-check for ana- 
lytical approximation schemes. 
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Figure 5: ++ spectrum beyond the critical Higgs mass; left: near to crossover, right: 
far from crossover. Triangles denote Higgs states, full symbols PF-ball states and circles 
Higgs states with admixture of pure gauge matter 
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